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We have studied structure formation in a confined block copolymer melt by means of dynamic 
. density functional theory (DDFT) . The confinement is two-dimensional, and the confined geometry 

f"^ ' is that of a cylindrical nanopore. Although the results of this study are general, our coarse-grained 

, molecular model is inspired by an experimental lamellae-forming PS-PBD diblock copolymer system 

(Shin et al, Science, 306, 76 (2004)), in which an exotic toroidal structure was observed upon 
' confinement in alumina nanopores. Our computational study shows that a zoo of exotic structures 

q_) , can be formed, although the majority, including the catenoid, helix and double helix that were 

also found in Monte Carlo (MC) nanopore studies, are metastable states. We introduce a general 
classification scheme and consider the role of kinetics and elongational pressure on stability and 
formation pathway of both equilibrium and metastable structures in detail. We find that helicity and 
three-fold connections mediate structural transitions on a larger scale. Moreover, by matching the 
remaining parameter in our mesoscopic method, the Flory-Huggins parameter x, to the experimental 
system, we obtain a structure that resembles the experimental toroidal structure in great detail. 
Here, the most important factor seems to be the roughness of the pore, i.e. small variations of the 
C/2 , pore radius on a scale that is larger than the characteristic size in the system. 
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I. INTRODUCTION 



Pattern formation of block copolymers in constraint situations or confinement is an important topic in polymer 
research, since the meso- or microstructure can be much better controlled when compared to the bulk. In general, 
O . molecular conformations and assembly are strongly influenced by confinement. In the absence of external constraints, 
the microstructure is dictated by the interaction between segments comprising the copolymer, the volume fraction of 
■ the blocks, and the molecular architecture. In a confined system, however, interfacial interactions, symmetry break- 
' ing, structural frustration and confinement-induced entropy loss play a determining role and may lead to structures 
that differ from the ones found in bulk. As a direct result phase separation in confinement has been the subject of 
extensive theoretical and experimental studies. [1] The intriguing prospects from a technological viewpoint are the 
novel structures that can be achieved and that may serve as scaffolds for other nanostructures. 
The simplest example of a block copolymer is a linear AB diblock copolymer. In the bulk several stable periodic 
microstructures can be formed, among which are lamellar, hexagonal and body-centered cubic phases. [2] The equi- 
librium behavior for AB diblock copolymers in the bulk has been mapped out both experimentally and theoretically; 
the situation for more complex multi-block and/or branched block copolymers is much less clear. Here, we focus 
on symmetric or nearly symmetric AB diblock copolymers. In the bulk, these copolymers microphase separate into 
lamellar microdomains with a characteristic equilibrium period L$, where grains of ordered lamellar microdomains 
k>( \ are randomly oriented. Global orientation of the microdomains can be induced by confinement of the block copolymer 
5—1 ' ^ n thin supported films or slits (one dimensional confinement). [3-6] Additionally, the interplay between surface fields 
(interaction with confining surfaces) and confinement effects (commensurability) can affect the phase behavior, and 
lead to the formation of surface reconstructions or hybrid structures. [7-9] 

A obvious next step in experimental and theoretical research is to consider systems where the confinement is effec- 
tively two-dimensional. This is the case when the melt is confined inside a cylindrical nanopore of radius R. From 
a conceptual point of view, this type of confinement differs from the ID confinement since: (1) there is only one 
confining surface, so there is only one surface field, (2) the nanopore solid surface is curved. By definition, surface 
reconstructions (analogs of the uni- or multilayer parallel lamellae Lu in slits) are therefore curved as well, which 
has an effect on the entropy contribution of the structure in the free energy. This can either result in a break-up of 
the nanopore in two relatively independent regions (close and away from the pore surface), or in curved multilayer 
structures that experience lamellar bending throughout the pore, with an increased bending towards the centre of 
the pore, (3) just like in slits, surface reconstructions may be able to adapt their layer spacing to some extent in 
order to suit the cylinder radius R. One can expect the influence of this type of frustration to be less than in slits, 
as there is no block-preference in the centre of the pore. However, for small R, structures can experience substantial 
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frustration due to incommensurability. To resolve this unfavorable situation, the lamellae can adapt an orientation 
perpendicular to the nanopore wall or the chains may find alternative packing and form hybrid structures, like in a 
slit, and (4) in addition to the commensurability issue for parallel structures or surface reconstructions (see (2)) the 
length of a nanopore can also affect the formation of perpendicular structures (incommensurability along the pore). 
In slits, these perpendicular structures are found for strong surface fields at incommensurable film thickness, and for 
weak surface fields independent of the film thickness. In nanopores, a mismatch of the pore length and the natural 
spacing Lo can lead to extensional forces, and may stabilise hybrid structures as well. 

Based on rather small data set of surface interactions and nanopore radii R, two pioneering computational studies 
(by Monte Carlo [12] and dynamic self-consistent field theory [13]) identified two structures: a slab and a multiwall 
tube morphology. [13] In analogy with the classification scheme for structures in a slit, [14,15] slabs can be associated 
with lamellae perpendicular to the confining surface (LjJ and multiwall tubes with surface reconstructions: lamellae 
that are parallel to the confining surface (Lh). One should note that this remark implies that the nature of these 
structures is the same as the lamellar structure found in bulk. It was concluded that two mechanisms control the 
structure formation in a nanopore : (1) in the case of weak surface interactions the lamellae orient perpendicular to 
the pore wall to form lamellar slabs and (2) for strong surface interactions, one of the blocks segregates preferentially 
to the pore wall, and the lamellae line up and bend to form concentric cylinders, the number of which is determined 
by the cylinder radius. The effect of commensurability was found to be less significant in this type of confinement, 
since no perpendicular structures were observed for incommensurate R and higher surface interactions. 
Very recently, the significance of these theoretical studies increased dramatically by the appearance of a number of in- 
tricate experimental studies of polystyrene- Wocfc-polybutadicnc (PS-&-PB) [16-19] and polystyrene- block-poly (methyl 
methacrylate) (PS-&-PMMA) [20] diblock copolymers confined in nanopores. In these studies, bulk lamella-, cylinder- 
and sphere- forming block copolymers were introduced into nanoscopic cylindrical pores in alumina membranes. For 
bulk lamella-forming block copolymers the predicted morphology [12,13] of slabs was not found, most probably due 
to the strong disbalance of surface-energies in the experimental setup (although these were not measured as such). 
The most frequent found structure was a concentric cylinders (multiwall tube [13]) morphology. However, also a new 
structure, a stacked-disc or toroidal-type [16] structure, was found inside a nanopore with an incommensurate pore 
diameter (d/L ~ 2.6, with d the pore diameter). 

As a result of these experimental findings, computational studies have considered the phase behaviour of confined 
symmetric diblock copolymers in more detail, by means of self-consistent field theory (SCFT) [21] and Monte Carlo 
(MC) [22-25] methods. Based on 2D calculations, the work of Li et al [21] aimed at constructing phase diagrams 
(X-W versus / = /a), and found lamellae for Ja = 1/2 (the diagram was for calculated for a single radius R — 8.5R g ). 
Although this work can be seen as a stepping stone for the understanding of the nature of phase transitions under 
varying conditions, their findings do not directly relate to the essentially 3D structures found in experiments. For 
instance, the slab morphology is out of the scope of their calculations. Moreover, the influence of the nanopore radius 
and the surface energetics was not considered in detail. Based on lattice MC simulations, Chen et al, Feng et al 
and Wang indeed found several new structures in 3D: a single helix, catenoid cylinder, gyroidal, stacked circle, and 
disordered structure in Ref. [22], a mesh, lamellae parallel to the pore axis, single helix and double helix structure in 
Ref. [23] and catenoid cylinder in Ref. [25]. A later study of Feng et al concentrated on the stability of the helical 
structure. [24] From a visual comparison we conclude that the mesh structure [23] is the same as catenoid cylinder, 
[22] . Similarly, one can claim that the gyroidal structure is in fact a defected structure, and comprises a coexistence of 
a mesh and a double helix structure. [22] We adapt a common notation, and conclude that this reduces the number of 
computed nanopore structures to seven: stacked disc, concentric cylinder, lamellae parallel to the pore axis, catenoid 
cylinder, disordered, single and double helix. 

Although new structure were found, a fundamental understanding of the underlying mechanisms such as achieved 
for thin films (effective ID confinement) is absent, except for the MC study of Wang [25]. However, the latter study 
is restricted to strong segregation, and most remarkably several of the other structures (for instance the interesting 
helix and double helix) were not found at all. Moreover, the origin of the experimental stacked-disc or toroidal-type 
[16] structure is still unexplained. More fundamental studies, including scans of a larger parameter space, are clearly 
needed for a deeper understanding. 

In principle, one can anticipate several regimes. Equilibrium morphologies are minima of the free energy, containing 
both energetic and entropic contributions. In the absence of any surface field and for a strong surface field, either 
one of these contributions is dominating, leading to the stacked disc and concentric cylindrical structures. Between 
these regions, for weak surface fields, the chains have more flexibility to adapt their packing and the system can 
and apparently will adapt other morphologies. Incommensurability may play a subtle role here. The details of this 
interplay remains to be determined, for instance by the calculation of a structure diagram in 3D, depending on R and 
the surface interaction strength. Since we consider lamella-forming systems, /a, the volume fraction of the A block, 
is fixed (in contrary to Ref [21]) and the diagram is conceptually similar to the one on Ref. [14] for confined films. 
Moreover, two factors have not been considered in detail yet: the influence of the kinetic pathway and the value of 
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XabN, where iV is the total length of the block copolymer. Although we will not focus into detail in the latter factor, 
from Ref [21] one can conclude that an increase of XabN may lead to a decrease of the number of concentric lamellae. 
Here, we focus on calculating the structural diagram and the question of stability. Whether interesting structures 
like double helix can actually be manufactured experimentally may be a very subtle issue. It is known that in strong 
confinement the dynamics may slow down and structures may be frozen into metastable states. This is particularly the 
case when the free energy difference between states is small or when there is a large energy barrier between different 
local minima of the free energy. In nanopores, this situation was actually observed in MC studies: in Ref. [22] for weak 
surface fields up to three different structures were found. Previously, a dynamic density functional theory (DDFT) 
study [13] identified the single helix structure as a long living intermediate structure between complete mixing and 
stacked discs for zero surface field. We will use this method here to consider the experimental system of Ref. [16] in 
detail. In contrast to traditional schemes of polymer phase separation dynamics where a Landau Hamiltonian is used 
with vertex functions calculated following the Random Phase Approximation (see e.g. [26]), we numerically calculate 
the free energy F of polymer system consisting of Gaussian chains in a mean field environment using a path integral 
formalism [27,28]. Our approach uses essentially the same free energy functional as in SCF calculations of equilibrium 
block copolymer morphologies by Matsen and Schick [29], but complements the static SCF calculations by providing 
a dynamical picture of the system. 



Here, we shortly discuss the DDFT method [27,30] for a bulk lamella-forming diblock copolymer melt. The diblock 
copolymers are modelled by a An a Bn b Gaussian chain (N = Na + Nb, fi = Nj/N). The confined geometry is 
a cylindrical pore with varying diameter R. Periodic boundary conditions apply in the x-direction, along the pore. 
Calculations are carried out on a cubic L X X Ly X L z grid with a spacing Ax that is related to the Gaussian bond 
length a via aAx^ 1 = 1.1543. [31] The spacing Ax is equal among different pore systems, making the free energy 
per volume element of Ax 3 easily comparable. Unless mentioned otherwise, all spacings are in units of this basic 
variable Ax. The pore is introduced into the simulation volume V by a masking technique; as a result the simulation 
volume contains both pore and mask points. [30] In mask points, elements of the subset V° = {r = (x,y, z) G 
V| |j (rc, j/, 2;) — [x , y{) , Z{))\\ > R} with yo = L y /2, zq = L z /2, all concentration and external potential field values are 
set to zero, except for the auxilary field Pm{ y ) = 1 that represents the mask itself. [30] The free energy for unconfincd 
systems is given by (see [27] for details) 



Here, k is the Boltzmann constant, T is the temperature, n is the number of polymer molecules in the volume V 
occupied by the system, and $ is the intra-molecular partition function for ideal polymer chains. The parameter k 
determines the compressibility of the system (the dimensionless «' = (inv = 20.65, with v the bead volume), and 
Pj is the mean concentration of the /-block (where the average is taken over V \ V ). For k — > oo the system 
becomes incompressible. The external potentials Ui and the concentration fields pi are related via the density 
functional. [28] The inter-chain interactions are incorporated via a mean-field with interaction strength controlled 
by the Flory-Huggins parameters xu ■ I n l me with our earlier work the interactions are specified by the parameters 
e°jj (in kJ/mol), [27,28] which are directly related to the Flory-Huggins parameters by xu = lOOOeJj/n^fcT (with 
ha Avogadro's number and T = 300 the temperature in Kelvin). In case of non-zero surface interactions an extra 
cohesive term is added to the free energy (1) equal to [30] 



In line with our earlier work, the interaction kernel is chosen Gaussian, and the important parameter e° IM denotes the 
scalar interaction strength [27] of bead I with the pore boundary. In the computations, e BM = indicating that the 



II. METHOD 





(1) 




(2) 



3 



S-blocks have no interaction with the wall, which is appropriate due to the fact that for an incompressible system of 
diblocks only the effective surface interaction £ = e AM — e BM is of importance. Using this we obtain 

F wa " = ]T e° IM F[ Pl , PM ] = &\pa, Pm] (3) 
i 

where T is 

F = cJ e-^ {r - r ' )2 p A {v)p M {v')dvdv' (4) 

v 2 

with C = (2^2)^ a normalisation constant and a the Gaussian bond length. Since the bead-bead interaction e AB is 
mostly considered a constant, we will use the notation (i?, £) to denote points in the structural diagram. 
For the simplest model, the evolution of the density fields is given by a Langevin equation [32] 

^M=MA M/ (r)+r?/(r) (5) 

with M a constant mobility, A the Laplace operator and r\ noise, distributed according to the fluctuation-dissipation 
theorem. Other transport coefficients with a more physical scaling behaviour, like the one for collective Rouse 
dynamics or reptation, exist, [32] but are in general too computationally demanding. We have recently found that 
our relatively simple model (5) can be appropriate to describe the experimental dynamics in detail. In Ref. [11] the 
experimental and calculated dynamics of cylinder and sphere forming diblock copolymers under an external electric 
field was compared, and good agreement was found based on DDFT with constant transport coefficients. In Ref. [10] 
the calculated dynamics of DDFT was shown to agree in full detail with SFM measurements of experimental dynamics 
in a thin film of a concentrated SBS solution. Finally, apart from the extra term in the free energy, confinement is 
accounted for by the boundary condition for the dynamic equations n ■ V/i/ = 0, with n the normal pointing into the 
solid object. [30] 



III. RESULTS AND DISCUSSION 

A. Nanopores: data from literature. 

Since the nanopore simulation data in literature is rather scattered, [12,13,22,23] we first present a short overview. A 
direct comparison of these different studies is complicated by the fact that important system parameters are different. 
A Monte Carlo (MC) method was used to simulate the phase behaviour for a diblock copolymer [23] with 

£ab = 0.3, 1.0 and 1.1 kT, and for a AiqBio diblock copolymer [12,22] with eab = 1-0 kT. In both studies all other 
interactions are zero: only the sas (in kT), the interaction between A-blocks and the surface, is varied. In the dynamic 
density functional theory (DDFT) calculations of Ref. [13], a AgBg system was considered for eab = 2.5 kJ/mol and 
varying e^s (in kJ/mol) (in that article, S is actually called M or mask). Using the formula in the method section 
for T = 300K , we obtain x = 1.0 and \N = 16 (weak segregation). Using the expression from Ref. [33], \ « 5(e/fcT), 
we can recalculate the parameters used in the MC studies: %iV « 15, 50 and 55 (weak to intermediate segregation) 
for Ref. [23] and \N ~ 100 in Ref. [12,22] (strong segregation), where it should be noted that Ref. [23] concentrates 
on \N ~ 15 (weak segregation). As a consequence, the results are distributed between weak, intermediate and strong 
segregation regimes. In Figure 1 we have combined the existing knowledge in schematic diagrams for different ^iV; 
we have only differentiated between stacked discs, concentric cylinders and alternative structures. It should be noted 
that the diagram for DDFT is essentially no phase diagrams, as the final structures are pathway dependent, and 
not minima of the free energy per se. This is a fundamental difference between static calculations aimed at deriving 
equilibrium morphologies, and our dynamic simulations aimed at mimicking experimental pathways, including visits 
to long-living metastable states. The important dimensionless spatial coordinate is the ratio of the pore diameter R 
and the lamellar domain distance L ; this parameter is on the vertical axis. On the horizontal axis is sas (in kT) 
of MC. Since the diblock is symmetric we can restrict ourselves to only positive values. We have used £,45 (in kT) 
w l/2e° AS (in kJ/mol) for the conversion of the DDFT to the MC value. 

We see that concentric cylindrical structures dominate all diagrams when £,45 is large, irrespective of the value of 
XN. In this case, the surface field dominates and one block is energetically favoured at the pore wall. For neutral 
pores, i.e. in the absence of a surface field (£,45 = 0), stacked discs dominate. In stacked discs the chains can adopt a 
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packing that is similar to the one in normal bulk lamellae. However, for the strongest segregated system stacked disc 
coexist with a single helix for this surface field. For small non-zero values of sas also other structures can be found. 
Although the data are sparse, one can conclude that structures with alternative packing dominate the diagram for 
large xN and relative small £^45. Incommensurability cannot be the only important factor since these alternative 
structures are even found for commensurate pore radii R/L = 1, [22] although there is some uncertainty since the 
bulk domain distance Lq was not determined explicitly in this study. We can only conclude from these data that 
especially weak surface fields in combination with strong segregation (high \N) lead to structures with alternative 
packing. 

B. Nanopores: system choice and boundary condition 

The system was chosen to model the experimental polystyrene- Wocfc-polybutadiene (PS-&-PBD) diblock copolymer 
in Ref. [16], which has a volume fraction of 0.56 for the butadiene block. The molecular model considered here is a 
S10-B12 (Ns = 10, Nb = 12) Gaussian chain. The interaction between S and B beads is chosen as esb — 2.0 kJ/mol 
or x = 0.8, and consequently xN = 17.6. A diblock copolymer melt of this molecular composition fs = 0.545 forms 
lamellae in bulk, and the bulk lamellar distance was determined as Lq — 8.6 (from here onwards, all distances are in 
units of Ax). Commensurability issues in the x direction, along the pore axis, may arise as a result of periodic boundary 
conditions. Since the effects of periodic boundary conditions have been considered in the literature, [35,36] we refer to 
these works for a detailed discussion. Here, we only note that the studies mentioned in the previous paragraph have 
not considered this effect in detail. Only for one dataset in Ref. [22], were both helices and stacked discs were found 
and no transition between the two different structures, it was shown that the number of MC simulations required to 
first find a helical structure peaks at particular cylinder lengths (our L x ). General conclusions from this study are 
not easy to make due to the absence of regularity. 

C. Nanopores: results and discussion 

We have considered pattern formation in a S10-B12 diblock copolymer melt confined in a nanopore of length 
L x = 32 = 3.72L , except for the larger pores where we have used L x = 34 = 3.95L and L x = 36 = 4.19i . 
In order to give a unifying description for the structural behaviour of slightly asymmetric diblock copolymers in 
nanopores, we first construct and discuss a diagram of simulated structures. Since these structures are the result 
of dynamic pathways and therefore may be metastable, we consider their stability by an interpolation procedure 
described below. Consequently, we focus on several factors that are important for structure formation: two types of 
incommensurability issues, due to packing frustration along the pore (associated with perpendicular structures) and 
perpendicular to the pore (associated with parallel structures), and the kinetic pathway, that may lead to arrested 
structures. In the second part, we discuss these issues in detail, and unravel the complex interplay of these factors 
and the surface field. 

Packing frustrations along the pore, originating from the requirement of periodicity at the two pore boundaries, have 
not received much attention in the past, except for the unconfined situation. [36] In this study, it was shown that 
this frustration, quantized by L x /Lq, determines the (in)stability of several perpendicular structure. This effect, that 
can be related to elongational stress or extensional force in an applied shear field, will be discussed in detail for our 
results. 

In the absence of extensional forces, the analogy to similar substrates in thin films suggests that the phase behavior 
in pores is due to an interplay of two factors: the strength of the surface field and confinement effects. [1,14,37] For 
vanishing or weak surface fields, the elastic chain deformation associated with parallel structures will be avoided, 
and perpendicular structures are favored instead. Due to the selective block-surface interaction, parallel structures 
(surface reconstructions) will be promoted for surface fields above a certain threshold field strength. However, the 
available space in the confined geometry dictates the degree of chain compression or extension required for parallel 
structures. Large chain frustrations give rise to an entropically unfavorable situation and may even prevent the for- 
mation of parallel structures. Alternative (in thin films, perpendicular) structures are then formed instead. In thin 
films, the surface fields were shown to be additive and affect only the first layers of structure. [37] Consequently the 
available space, quantized by R/L in pores, is important for both factors. The surface field and frustration due to 
confinement will strongly influence structure formation in small pores. In larger pores, the surface field has a limited 
range, and possible structure frustration can be distributed over more layers of structure. 

We note that the diblock copolymer considered here is slightly asymmetric. As a consequence, the structure diagram 
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is not completely symmetric as well. In the absence of surface interactions, the shortest (S) part of the chain is prefer- 
entially found close to the surface due to entropic effects. The value for which energetic and entropic contributions are 
balanced shifts to a small positive £. [14] Moreover, for parallel structures in the curved geometry, this chain asymme- 
try leads to curvature effects, since the spontaneous curvature associated with each of the different blocks is slighlty 
different, and asymmetric packing frustrations due to the different S and B domain sizes. In contrast to thin films, 
the majority component of the centre layer in a parallel structure is not prescribed by surface energetics. Packing 
frustrations are therefore expected to be less significant than in thin films, where in most cases the block next to the 
surfaces is prescribed by the surface fields. However, the number of domains/layers for each of the blocks can differ. 
To focus on this effect we adapt two notations: if the majority component in a parallel structure \sS — B — S — B — S 
or B — 5 — B — S — B, on a line perpendicular to the pore axis from wall to wall, we denote the structure as L| 1 ± (for a 
concentric cylindrical structure, in line with the slit notation [14]), and by the order of the majority component in the 
layers from wall to centre, SBS or BSB. Note that here the total number of S and B domains differ (2 or 3), as well 
as the composition of the central cylinder, associated with the largest curvature. Curvature is an important factor and 
a complicating factor, when compared to the thin film situation. In principle, each of the layers in a parallel structure 
can adapt their thickness to some extent in response to global mismatches. Simple volumetric arguments indicate 
that the domain spacing in nanopore confinement depends on the local curvature, and therefore on the absolute radial 
position of S — B interfaces in the pore. [25] This is particular the case in the limit of strong segregation; for weaker 
segregation, the situation may be different, since the blocks are somewhat miscible. As a result, the spacing may 
heterogeneously deviate from the bulk domain spacing (Lq = 8.6, so D$ — 3.9 and Db = 4.7), near the pore centre, 
where the curvature is highest, and in the wetting layer, where confinement effects are most severe. For a symmetric 
diblock copolymer and larger radii R, we previously showed that packing frustrations in parallel structures are mainly 
relieved by rearrangement in the cylindrical centre region; in the layered structure away from the centre, the chain 
conformations are rather unaffected by the local curvature. [13] We will consider this interplay of surface field and 
frustration in detail for our asymmetric block copolymer. 

We have varied the surface interaction in a range of negative and positive values £, between vanishing (non-selective) 
to intermediate (selective, S or B) surface fields. A range of small R up to approximately Lq was chosen to consider 
the details of this interplay in strong confinement. Moreover, also a few larger R were considered to study these 
effects separately. The structures are displayed in Figures 2-4, and the free energies associated with these structures 
in Figure 5. In the selected region in the (R, £) diagram, we expect that confinement effects, surface fields and/or the 
interplay between these two gives rise to alternative structures. In the remainder, we use the relevant dimcnsionlcss 
parameter R/L instead of the bare radius R. For all calculated final structures, the free energy was monitored and 
remained constant. 

We introduce a short-hand notation for the perfect structures: catenoid cylinder are denoted by PL or PL{I) (de- 
pending on the sign of £, the majority component I of the structure is either S (negative) or B (positive)), disordered 
by D, single helix by H and double helix by DH. The lamellae parallel to the cylinder axis are perpendicular to the 
pore wall but also parallel, and we denote them as In some cases, defected structures are remarkably stable. 

Based on visual inspection, and knowledge about metastable intermediate states along the pathway of formation of 
the perfect structures, we assign a symmetry group preceded by the letter 'c£' (for instance, dDH is a defected double 
helix). Other new structures will be discussed and annotated as they appear. 

Focussing on the kinetic pathway, we remark that all structures in the diagram of Figure 2 were obtained following 
a diffusive pathway (equation (5)). Earlier work showed that microphase separation in confined systems often starts 
close to the pore wall. [38] For the pore dimensions considered here, this effect will be small, but, especially for 
stronger surface fields, surface reconstructions will initially form, resulting in overall coverage of the pore wall by a 
single component. The transition to more stable structures with a different surface coverage (for instance, helices or 
stacked circles) requires considerable transport of material away for the wall, into the center of the nanopore. The 
fundamental mechanism of this transition is important, as it may hint when and where this process may be kinetically 
trapped. First, we consider which of the calculated structures is an equilibrium structure. For varying £, the free 
energy (see Figure 5) associated with a particular structure (defined by ps', due to incompressibility, ps is then also 
fixed) is given by (see also (1)) F tot = F + F wal1 = F[p s ] + £F[ps,Pm\- Since both F[p s ] = a and T[p s ,p M ] = b 
have a constant value for a particular ps (pm is the mask field, and fixed by definition) the free energy values for this 
structure and varying surface field can be found on a line defined by a + £b. We only need two datapoints (or the 
values of a and b) for the same structure to determine this line. However, we have to be careful using this procedure 
as structures with the same morphology type may differ in detail (and therefore in the values of a and b), such as 
geometrical quantities and the degree of segregation. For instances, interpolation suggests that the L± structures for 
£ = —0.4, —0.3 and £ = 0.1, 0.2 (both for R = 0.70L ) arc different, although there is no structural difference in terms 
of easy computable geometrical quantities. Only when we consider the bare density values wc find that the maximum 
of the concentration field #s( r ) f° r £ = —0.4, —0.3 is slightly higher than for the equivalent structures at £ = 0.1,0.2. 
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1. Structure diagram and stability 



The structure diagram is shown in Figure 2, and is indeed not symmetric. Equilibrium structures are distinguished 
by a grey background; this determination is based on the interpolation approach described above. Considering the 
general features of this diagram, one observes that several structures in this range are metastable. For large R/L , 
only Li and L±_ are stable; for small R/Lq also other structures can be stable. The effect of incommensurability is 
limited, since we do not observe outliers of perpendicular structures for the larger radii considered (where the effect 
of the surface field is relatively small). However, for positive £ (pore surface likes B) and R/Lq = 1.28 or 2.21, the 
i|| structure is only metastable and the transition to equilibrium structures trapped. Overall, this suggests that the 
pore radii R/Lo = 1.28 and 2.21 are incommensurate, while R/Lq = 1.74 is commensurate. Moreover, we see that 
for negative £ (pore surface likes S), Lu is stable for £ < —0.4, independent of the radius R/Lo. This asymmetry is 
due to the curvature effects and asymmetric packing frustration mentioned above. 

Focussing on specific structures, we see that for large absolute values of £ surface reconstructions, PL, dPL, dLo, 
and Lu, dominate. In particular, we find PL and (SB or BS) for R/L G {0.58,0.70}, Limi (SBS or 
BSB) for R/L e {0.81,0.93,1.05,1.28}, L\ l2 (SBSB or BSBS) for R/L = 1.74, and L [[ 2 x (SBSBS or BSBSB) 
for R/Lo = 2.21. An analysis of domain spacing is given in Table I; the values are derived from the interface locations. 
Disconnected concentric cylinders dL\\, where the cylinder close to the pore surface is broken up into 4 disconnected 
parallel lamellar patches, only appear for R/Lo = 0.93. The distance between the patches increases with increasing £ 
in this region. The Lj_^ is the same as the structure in Fig 8b in Feng et al, [23] and related to dL^ (see discussion 
lateron). The number of perforations in PL can vary and show hexagonal ordering when considered in the 2D plane. 
Geometrically, the PL structure can be seen as an intermediate between L|| and H. Defected perforated lamellae 
dPL, distinguished from PL since the perforations do not show ordering on a larger scale, are found on the boundary 
between parallel and perpendicular structures. In response to the perforations the central cylinder sometimes adapts 
an oval cross-section. Defected PL only appear for negative £ (surface preference to S). The majority of surface 
reconstructions (PL or L\\) for positive £ (surface preference to B) are not stable structures. They are examples of 
kinetically trapped structures, due to the presence of the pore surface, that gives rise to L\\ related structures in the 
early stages of phase separation. 

Perpendicular structures, titled tacked discs L^'* and L± (4 discs for all L x considered) , are found in the centre of the 
diagram, for small £. For small relative radii R/Lo, the discs are perpendicular to the pore surface, L±. For larger 
R/Lo there are only three discs, and they are tilted with respect to this surface, L^'*. For the largest pore radii, 
R/Lq = 1.74 and 2.21, the pore length L x (L x = 34 = 3.95L and L x = 36 = 4.19-Lo, respectively) is somewhat larger 
than the standard L x = 32, and the structure is Lj_. The L^/L^ structure are two coexisting perpendicular L*^* 
structures, with opposite tilt angles. The details of these mixed structures, not representing equilibrium structures 
for obvious reasons, will be considered later. 

Helical structures, H , DH and a new GH structure, are found at the rather broad boundary between perpendicular 
and parallel structures. These structures possess features relating them to both perpendicular and parallel structures: 
the axis of winding is parallel to the pore axis, and the wall-coverage is nonuniform. The GH structure in (1.05, 0.4) 
is related to the helical structures, but differs topologically since small helical patches are three-fold connected with 
helical patches on the opposite side, and these connections join into two cylinders parallel to the pore axis. The GH 
for (1.16,0.5) and (1.40,0.4) share this property, but instead the connections form two lamellar patches and three 
cylinders, respectively. We call this structure gyroid- helical (GH) because of the three-fold connectivity. Stable helices 
(left- and righthand) and a double helix are formed for R/L < 1 and negative £. The helical structures for positive £ 
are all metastable. Metastable GH are only found for larger effective radii, R/Lq > 1, adjacent to both parallel and 
perpendicular structures. 

Several coexisting structures, Lj_/PL, L^/H, L*^* /GH, L^/H and L±/L\\, and defected structures, dL 1 ^ 1 , dL^, 
dGH , dLj_ ii are found directly adjacent to their 'perfect' counterparts, and stay defected after many TMS. Apparently, 
the driving force for the removal of different types of defects in these structures is rather small, and as a result the 
structures are kinetically trapped. The only exceptions, coexisting H/PL and PL/Lu for positive £ that are not 
adjacent to a PL structure, show that PL is associated with a free energy close to the one for H and Lu, respectively. 
For R/Lo = 0.81 and —0.4 < £ < 0.0 the structures are very defected, and remain as such, even after a large number 
of extra timesteps (TMS). 



2. The value of L x /Lo: elongational stress and perpendicular structures 

Earlier work for a 'soft' confined system [36] concluded that the L±, L 1 // 11 and H structures are related. Their 
stability depends on the extensional force on the system, originating either from an external field (experiments) or 
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from boundary conditions (computations). Extensional forces may be present in experimental and computational 
studies dealing with structure formation in nanopores, except for a computational study in the strong segregation 
limit (SSL), [25] where MC based on a grand canonical ensemble was employed, but have not received much attention. 
Especially in computational methods considering a canonical ensemble and a single L x , this type of commensurability 
issues cannot always be avoided. We argue that a deeper understanding of this effect is relevant. In experiments, for 
instance, extensional forces play a role when the pore surface is very rough, the pore length is very small (often the 
case in applications considered in soft nanotechnology) , and when shear fields are present, as is the case in a recently 
developed experimental technique for the fabrication of nanowires. [39] In the present study, we have considered this 
effect by varying L x for a small set of selected parameters: L x — 34 = 3.95Lo for (0.58,-0.1) and (1.28,0.1) and 
L x = 36 = 4.19i for (0.70, -0.5), (0.81,0.2), (0.93,0.2) and (1.05,0.1). In all cases the microstructure evolves into 
Lj_ (4 discs), with an associated free energy that is lower than for the orginal structure for L x — 32. A detailed 
analysis of the structure evolution revealed that L*^* and H structures can also be kinetically related to L± : in some 
cases, they mediate large-scale structural reorganisation. This phenomenon can be observed from the helicity that 
appears during the evolution of the L± structure in a pore of L x = 36 and (R/Lo,£) — (0.70, —0.5) (see figure 8 for 
the details of the formation pathway for different L x ). It is in agreement with earlier findings for a fully symmetric 
AB diblock in Rcf. [13]. 

Apart from the single additional calculation for each of pores and L x ^ 32, we used the L± structures to compute the 
free energy for an additional surface field strength £. Interpolation between these two values enables us to reconsider 
the stability of the structures in Figure 2 with respect to Lj_. We make two general remarks: i) Although the 
interpolation technique is valuable for comparing the stability of different microstructures in pores of equal length (for 
varying £), in principle this procedure cannot be used to differentiate between microstructures in pores of different 
length. For instance, an instantaneous change of the pore length (or equal, extensional force) during the evolution 
gives rise to a deformation of the structure, and possibly a transition to a more stable structure. However, in contrast 
to self-consistent field (SCF) techniques, DDFT does not impose symmetry, and changes in free energy associated 
with a deformation of an existing structure, for instance an affinc deformation of a helical structure, cannot be 
quantized directly. Instead, additional calculations are necessary for many alternative structure types and are very 
time consuming in general. [36] We therefore anticipate on the results of the SSL study, [25] that considered nanopore 
structure formation in the absence of clongational stress. This study identified, besides lamellar structures, only 
stable PL. Consequently, we adapt an practical approach and use all data obtained by simulation (see Figure 5), 
independent of L x , to suggest a phase diagram of equilibrium structures (see Figure 6), without considering the 
deformability of the non-lamellar structures. We note that the stability of the PL phase could not be determined 
indefinitely, since the alternative L| phase was not formed for these small radii. Only for R/L = 0.58 and positive £, 
concentric cylinders (BS) were formed, and the region where PL is equilibrium structure could be determined, ii) The 
additional simulations for L x =/= 32 were carried out for relatively weak surface fields, where perpendicular structures 
are likely to form. Interpolation indicates (see figure 6) that the equilibrium structure can be perpendicular even for 
relatively strong surface fields, for instance due to incommensurability along the radial direction. However, the large 
structural rearrangements required for the transition to (stable) perpendicular structures, starting from the (unstable) 
parallel structures that are initially formed due to the strong surface interaction, may lead to arrested structures, i.e. 
microstructures that are kinetically trapped in metastable states. A good example is the DH structure in (0.81, —0.5). 
Although the phase diagram of Figure 6 shows that Lj_ is the equilibrium structure, we challenged the stability of 
the DH structure by choosing a number of systems with different pore length L x (each system was quenched from 
a homogeneous mixture at TMS=0, see Figure 9-10). We never obtained the equilibrium structure as a result of 
the dynamic pathway. Instead, upon a variation of the extensional pressure, we obtain structures with an overall 
parallel orientation, that lack symmetry on a larger scale. Only for L x = 32, 35 and 38 we find an almost perfect DH 
structure. We rank these other structures in three different classes, labelled as DHi. Structure DH1 can be seen as a 
highly interconnected and defected structure, where the centre lamellae contains large holes and is connected to 
the other structure. Structure DH2 is completely symmetric; the axis of symmetry is L x /2. Structure DHi consist 
of two disconnected sheets of different structure, one similar to PL(B) and the other similar to PL{S). For L x = 64 
the structure is very defected and a combination of other structures. Comparing the free energies in Figure 10, we 
find that the one of the metastable DH structure (L x = 32) is the lowest, indicating that the extensional force is 
minimal for L x = 32 and the structure is trapped due to the diffusive kinetic pathway. The finding of DH for other 
L x shows that the factor 3/L (with L the particular spacing of the structure) plays an important role. In general, we 
see that an increase of the nanopore length slows down the separation dynamics. Close examination of the structures 
in diagram 2 reveals that the calculated structures for —0.4 < £ < 0.0 are defected DH or resemble the intermediate 
DHI. In all cases, the kinetics disables the mass transport necessary for the DH — > L± transition. 
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3. The value of R/Lq: commensurability and parallel structures 

From the phase diagram in Figure 6 (see discussion for the derivation of this diagram above) we distinguish two 
specific values: R = O.8IL0 and R = 1.7ALq. For R = 0.81Loj the pore radius R < Lo and the surface field is therefore 
strong, but confinement effects apparently prevent the formation of parallel structures, in favor of the L± structure, 
up to large surface field stengths. Since this radius also marks the transition of SB (0.70i ) to SBS (or BS to 
BSB), we conclude that R = O.8IL0 is incommensurate. For R = 1.74Lo > L$, the surface field influence is much 
less. Apart from a symmetric Lm region for large surface fields (we note that the number of S and B domains is the 
same, independent of the wetting block), only a small region of perpendicular structures is observed for almost neutral 
pores, and we conclude that the pore size is commensurate. We consider the simulation results in Figure 2 in more 
detail. Domain distances are shown in Table I, and are inexact for the PL structure, which is not radial symmetric. 
For the smallest R = 0.58Ln 7 parallel structures are found for strong surface fields. Most of these structures are 
PL (PL(S) or PL(B)), and L\\ is only found for (0.58,0.7). The interpolation procedure showed that, at the onset 
of the L|| region, the free energies associated with both parallel phases are relatively close (the value of b is almost 
equal, see discussion of the interpolation procedure before), which explains that Lu is found for £ = 0.7 due to kinetic 
trapping. Upon comparing the PL and L\\ structures geometrically, we see that the formation of necks on the central 
cylinder in PL leads to an increased curvature of the S — B interface and contact area of the non-preferred block with 
the pore surface. If we compare the position of this S — B interface for this incommensurate and the commensurate 
situation (R/Lq = 1.74, S — B interface closest to the pore centre, see Table I) we observe a reduced thickness for 
positive £ and equal thickness for negative £. This illustrates the curvature effect mentioned before, meaning that 
the elastic chain deformation associated with the formation of can much easier be facilitated when the interface is 
curved towards the shortest S'-part of the chain, than to longest B-part. The same phenomenon can be observed for 
R = 0.70Lo. Here, concentric cylinders {L\\) are completely absent in the simulated range. Alternative PL structures, 
very similar to the ones for R = 0.58i apart from the number of perforations, can be found, but equilibrium PL(S) 
structures are only found for rather low £ values, again due to the curvature effect. In order to fill the pore, the 
domains are rather extended (Table I). The formation of parallel structures requires strong stretching of the part 
of the chain that contains the wetting blocks, giving rise to an entropic penalty that is only counterbalanced for 
stronger surface repulsion. For R = O.8IX0 surface reconstructions L^ A i (SBS or BSB) are found and PL is absent, 
apart from the coexisting H/ PL structure (£ = 0.5). However, from the phase diagram, we conclude that all surface 
reconstructions L\\ are metastable. and formed due to kinetic factors. The analysis of the domain distances in Table 
I shows rather small domain sizes, and therefore strong compression of the chains. Although L^i is only found 
for positive £ (R = 0.93-Ln)j we observe many other parallel structures, suggesting that this radius R is somewhat 
commensurate. However, from the analysis in Table I we see that the domain spacing for this radius is still rather 
small (compared to bulk spacing, and the spacing of D3 + D4 = 9.0 found for the inner two layers in the commensurate 
BSBS structure for £ = 0.7 and R/Lq = 1.74), especially towards the center. Morover, the equilibrium structure in 
this range is L± (see Figure 6). We focus on the disconnected lamellar patches (dLn) for £ < —0.2. Because of the 
masking procedure, the cross section of the pore is not completely circular, and the curvature of the pore wall varies. 
Hence, one could expect that metastable concentric cylinders and formed due to kinetic factors, and that the origin 
of these disconnections, appearing in regions of high curvature, is the chain-stretching required for the formation of 
undefected Lm. Close examination of the kinetic pathway towards the Lmi structure (figure 7d) shows dLu as an 
intermediate structure. In particular, along the pathway two of the lamellar patches and the central cylinder have 
merged to form the centre lamellae. In the early stage, undulations lead to connections between the different structural 
elements in dL\i (very fast). Subsequently, these connections merge (and others are removed) by a complex process 
that includes the formation of holes (very slow). The structures along this pathway resemble the metastable DH1 
structures that are found in (0.81,-0.5) for varying L x . Finally the Ljii structure forms. Apparently the removal 
of defects depends on the surface field strength and is slowed down considerably by an increase of £ from —0.1 to 
0.0. In MC calculations [23] the L±^ structure was also found for a pore radius R w L . It was suggested that it 
represents an intermediate between the L|| and L± structures, but no further explanation was given. The stability 
analysis suggests that all non-perpendicular structures are metastable. We claim that this structure represents a 
kinetically arrested alternative for stable L±. Both structures are formed for small to vanishing surface fields, in order 
to avoid the elastic chain deformation due to the curved S — B interfaces in concentric cylinders. Following the same 
reasoning, the dL\\ structures (instead of L||) for £ < —0.1 are metastable structures along the pathway to this Lj_,|| 
structure, and kinetically trapped since the undulations necessary for the formation of this structure are surpressed 
by the strong surface field. 

For R — 1.05L the situation is weakly incommensurate, and L± is the equilibrium structure between —0.4 < £ < 0.8. 
This can also be seen from the slightly perturbed spacing in Tabic I for £ = —1. For —0.5 < £ < —0.3 in Figure 2 
(dPL) holes are formed in the region of highest curvature, and we find that the density field in the centre of the pore 
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is slightly higher for £ = —0.5 than for £ — — 1 (no holes). This dPL structure was also observed as the initial stage 
in the formation of DH 1 H and L± structures for these £ and smaller R, and for the slow transition to dLff* for 
£ = —0.2, —0.1, where we simulated up to 40000 TMS. For positive £, the parallel structures are also metastable, but 
the free energy associated with the PL(B)/Lu coexisting structure (£ = 0.6) is slightly below the one extrapolated 
for perfect Ln. This is once more an indication that the chain deformation associated with concentric cylinders can 
be relieved by the formation of necks, and therefore PL is more stable. Moverover, this transition from L\\ to PL 
(and eventually the transition, via other intermediate structures, to L±) is obviously very slow. The parallel structure 
Lj_ M is absent, showing that the kinetic pathway leading to this metastable structure is rather unique. 
For R — 1.28.Lo, the phase diagram suggests that the pore size is slightly too large for Lpi . The detailed analysis 
in Table I shows that the chains in the center of the pore are stretched, compared to the bulk domain thicknesses 
(Db = 4.7 and D$ = 3.9), and to the L^ ^i structures for other R. The boundary for the equilibrium L± phase shifts 
to stronger surface fields, but this shift is small compared to R — 1.05Lo. The larger shift for positive £ is due to the 
curvature of the S — B interface of the central cylinder, which is curved towards the shortest block for negative £. 
For R = 1.74L surface reconstructions Lu 2 are found for positive and negative £. Moreover, only at the boundary 
with defected structures, these surface reconstructions are metastable. The stability region of Ln is almost symmetric 
around £ = 0.1, and we conclude that the situation is commensurate. The chains are slightly stretched (see Table I) 
compared the spacing to the domain spacing in the bulk due to curvature effects. A comparison of the spacings in 
pores for R/Lq = 1.74 and 2.21 suggests that especially the chains in the centre of the pore are compressed for the 
largest pore. The driving force for the transition of the mestastable surface reconstructions to L± is rather small (see 
Figure 5). 

4- Detailed kinetic pathways 

Although some kinetic pathways have been considered in the previous subsections, we focus on a few more repre- 
sentative pathways here. For R = 0.70i , different kinetic pathways were found leading to the final L± structures 
for £ = —0.4 (figure 7a) and £ = 0.1 (figure 7b). For the negative £ = —0.4 the initial structure is L||. At later 
stages (200 TMS) the cylinder in the centre of the pore has broken up to form short and disconnected cylinders, while 
part of the outer curved lamellae stays intact. Remarkably, parts of this outer curved lamellae remain intact during 
the transition and mediate the mass transfer from the wall into the pores interior. This is the source of the toroidal 
structures that appear as intermediate, at the positions where the stacked discs are later formed. For the positive 
£ = 0.1 the initial structure is bicontinuous (or very defected PL(S)) with many three-fold connections in the centre 
of the pore. These connections can be seen as defects, as very early in the evolution there is an excess of material 
at positions where stacked discs are later formed. Four-fold connections (small lamellar patches) and helical domains 
mediate the mass transfer needed for the formation of perfect stacked discs in this case. 

The most interesting structure for R = 0.81L is the double helix (DH) for £ = —0.5. Figure 7c shows that the 
formation of this structure is a rather slow process, and starts from L\\ (100 TMS). In the early stages, equidistant 
perforations form in the outer curved lamellae, similar to the ones formed for R/L = 1.05 and —0.4 < £ < —0.3. 
These perforations grow in time while the cylinder in the centre of the pore undulates and breaks up to form a 
bicontinuous network. In contrast to before this bicontinuous network itself forms the nucleus of the final structure, 
which is the double helix structure. Due to its helical nature, defect removal cannot be mediated by the formation 
and breakage of helical connections and thus the process is rather slow. 

The metastable parallel structure GH is only formed for positive £, and either stabilized by an extensional force or 
kinctically trapped along the pathway of the Lh — » L±_ transition. Its finding adjacent to perpendicular structures 
(deformed L±) and Lu for R = 1.28Lo ; or perpendicular structures and PL, coexisting with Lu and H, for smaller 
pores suggests that it may evolve from PL. Close examination of the pathway for R/Lq = 1.28 (£ = 0.4) and 
R/L = 0.93 (£ = 0.3) shows that the initial structure L\\ rapidly transforms (after 100 TMS) into a very defected 
structure, where the perforations in the outer B layer are uneven and randomly positioned. Since PL only appears in 
coexistence with other structures for positive £, we conclude that the kinetic pathways for the formation of GH and 
PL must be different. 

For R = 1.28£o, we concentrate on the effect of elongational stress on the evolution of the L±/L± structure, 
where the equilibrium structure is L±. In particular, we consider the structure evolution for £ = 0.2 in detail, since 
incommensurability issues with respect to the pore radius are absent. First, we note that the cylinder in the centre of 
the pore for £ = —0.2 show some undulations, but the number of pores in the outer layer remain constant for longer 
simulation times, and no connections with the inner structures are formed within the limits of simulation time. These 
connections are required for the transition/nucleation of one phase (in)to the other. For £ = —0.1 a defected structure 
without apparent symmetry at TMS=10000 evolves into a helical phase at TMS=25000. Since the axis of winding 
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is away from the centre, we conclude that the structure can be seen as a majority under some angle (with the 
pore wall), connected to a minority L^'* under a reverse angle. As a natural consequence, the connections between 
these two partial structures have a helical nature. We observe similar structures for £ = 0.1 and £ = 0.3. For £ = 0.2 
the axis of winding coincides with the centre of the pore, and we could classify this structure as analogous to the DH 
structure for smaller R. However, in line with the other structures for this radius this could be seen as a merger of 
two equally sized clusters of oppositely tilted and connected stacked disc morphologies. Moreover, in the evolution for 
£ = 0.2 (see figure 7e) we observe a fundamental mechanism: the orientation with respect to the pore wall of the total 
structure, as well as the change of chirality of local structure, are mediated by melting and reconnection of helical 
sections. In this way, unfavorable orientation and winding in part of the structure (for instance, giving rise to holes 
or increased curvature) can be avoided. We conclude that all these structures —0.1 < £ < 0.3 are actually coexisting 
jjilt unc [ er two different angles. The helical connections between these two cluster structures are relatively stable 
in time, and analogous of the double periodic array of saddle surfaces found in lamellar forming thin films under an 
electric field. [40] Apparently, they do not cost much energy to the system. This is confirmed by the free energy plot 
in Figure 5, were the free energies associated with these different helical structures fall almost on a single line. 
To see whether the elongational stress is the cause for this fundamental mechanism, we also considered the structural 
evolution in the absence of incommensurability along the pore, £ = 0.3 and R = 174L (see figure 7f). The detailed 
analysis confirms that the fundamental mechanism identified for formation of structure under elongational stress, i.e. 
reorientation on a higher structural level by helical connections, is also instrumental for the formation of equilibrium 
L± structures. One can clearly see that the tilted and interconnected discs that form in the initial stages (1000 
TMS) slowly transform into perfectly stacked discs or L± (25000 TMS). The mechanism is a moving front of helical 
connections, just like shown in Figure 7e. This effect can also be observed for R = 2.21L an d £ = 0.0 or 0.2 (not 
shown), where the initial coexisting structure Lf 1 / Lf * (TMS=10000) evolves into L± (TMS=20000). From these 
observation we conclude that coexisting stacked tilted discs with helical domains are the preferred transient structure, 
independent of the elongational stress. The most probable cause is that interconnected tilted structures can much 
easier accommodate the necessary large scale structural rearrangements that are necessary for the formation of L±. 
Helical connections play an important role in the removal of defects, and slowly change the tilt of the intermediate 
structure, leading to the formation of equilibrium L± structures. 

D. Discussion: comparison to other studies 

1. Experiments 

Up to now, experimental groups have considered structure formation of (nearly) symmetric polystyrene- block- 
polybutadiene (PS-6-PB) [16-19] and polystyrene-Wocfc-poly(methyl methacrylate) (PS-6-PMMA) [20] confined in 
nanopores. Almost all observed experimental microstructures belong to the class of concentric cylinders. This does 
not come as a surprise, as for most diblock copolymers one of the blocks has a higher affinity for aluminum. We have 
previously concluded that, especially for the larger pore sizes considered, surface fields above a certain (relatively 
low) threshold value lead to so-called surface reconstructions: concentric cylinders. In the experimental procedure, 
the block copolymers were introduced into the alumina nanopores at elevated temperatures via capillary action and 
subsequently cooled. After this procedure a weak base was used to dissolve the alumina and produced the free 
standing rods of block copolymer that were imaged. The thus obtained rods show several features: the cross-sections 
vary in shape (from ellipsoidal to spherical) and the thickness of the rod varies substantially along the rods (due to 
the roughness present in the alumina pores). As a result, the actual pore radius can only determined approximately. 
Although some measurements of surface tension were performed, exact numbers quantifying the affinity of PS, PB 
and PMMA for the aluminium pore material are missing. 

We conclude that the stacked-disc or toroidal-type structure in nearly symmetric polystyrene- Wocfc-polybutadiene (PS- 
fe-PB) block copolymer [16] is the best candidate for a true comparison between our calculations and experiments. 
We therefore choose the simulation parameters as close as possible to the values for this experimental system. The 
alternative structure was found at a temperature of T = 398. 15K, above the glass transition temperature [17,19], in 
a system with a molecular weight of M — 18400 (the number of monomers N ~ 250). The experimentally measured 
repeat period Lq = 17.6nm, giving rise to an incommensurate pore radius R/Lq ~ 1.3. We note that the effective 
FH parameter \ depends to some extend on the details of the experimental system. Using earlier measurements by 
the same group we find \ — 0.058 [41] and obtain \N « 14.5. For a coarse-grained chain with N = 22 (A10-B12) 
we therefore find an effective Flory-Huggins parameter xsb — 0.66. Here we consider xsb = 0.6 instead, and focus 
on the influence of reduced FH parameter in general (the calculations in the previous sections were carried out for 
Xsb = 0.8). Although we use the bare FH parameters in the calculations, the referring eg B = 1.5 kJ/mol for the 



11 



choice of the temperature (T = 300K) in the standard input in DDFT. 

The results of calculations for xsb — 0.6 and two radii R (and varying £) are shown In figure lla-b. We conclude 
that for the largest R = 19 the effect of reduced \ is rather small. Comparing the results for R = 19 in figure 11 
and R = 2.21L m figure 4 (for the same cylinder length L x ) we see that the phase boundaries do not been shift 
much. For low and high values of £ concentric cylinders are found with an additional layer when compared to the 
results for higher block-block interaction. This finding is well in agreement with the 2D SCF calculations of Li et al 
[21] where a transition from L$ to L2 (notation of Ref. [21]) is found for x^V > 13. The tilt of the L± for small £ 
reflects the different bulk lamellar distance L associated with the change in xsb and is due to an extensional force. 
For the smallest R — 9, concentric cylinders {Lu) are found for £ < —0.1 and £ > 0.5, with the same number of 
layers as for larger xsb- For 0.1 < £ < 0.3, a perfect L± structure is obtained. Constructing the lines associated with 
different structures shows that £ = 0.0 and £ = 0.4 are indeed located on intersections of different lines. Based on 
the experimental knowledge that the B block is preferentially found close to the pore surface, we further concentrate 
on £ = 0.4. For R = 19 the structure contains perforations, which will not be visible in the experimental imaging 
procedure. This structure will appear as concentric cylinders (L\\). For R = 9, the structure is a mixture of Lu and 
L±. The standard element are two opposite small lamellar patches oriented parallel to the pore surface. Along the 
pore, only the orientation of the elements changes and neighbouring elements are orthogonal (see figure 12a). The 
fourfold connections (instead of three) between the elements are short and curved cylinders. This mixed structure 
has interesting features that somewhat resemble the experimental toroidal-type structure of Shin et al [16]. However, 
the structures still do not match in all details. Additionally, we varied the small pore radius for fixed £ = 0.4 (figure 
11c). This leads to Lu for R — 5,8 and L± for R = 6, 7, 10 and 11. In general, we conclude that the decrease of %iV 
reliefs some of the stresses present in the stronger phase segregated system, and leads to faster relaxation towards a 
stable structure. 

Finally, the small scale of the experimental nanopores leads to constraints that are not effective in other systems. 
The results of the previous section show that the structure- type can be very sensitive to the pore radius. Small 
variations of the pore radius (due to roughness) along the pore at relatively larger length scales may therefore affect 
the structure formation to some extent. This is in particular the case in small nanopores, where slightly different pore 
radii will promote different structure-types and connectivity issues (of different types and periodicity) may appear. 
To study this effect, we prepared a pore with two pore radii, R = 9 and 10 (£ = 0.4), equally distributed along the 
pore (the thickest part is located in the middle of the pore). Using L = 8.6, the bulk characteristic lamellar spacing 
determined for the system with higher xsb, the relative pore sizes R/L vary between 1.05 and 1.16. Keeping in 
mind that Lq will decrease for decreasing xsb, this range matches the experimental value of 1.3 reasonably well. 
Following the evolution of this structure (see figure 12b), we initially observe that simply the stable structures for the 
individual pore radii are formed and connect (Lu for R = 10 and the cage structure for R = 9). However, eventually 
the interface between the different structures acts as a driving force 1 for a transition to another structure. The final 
and stable structure for the 5* block is completely perpendicular, with stacked disc for R = 10, but further away from 
this region the stacked discs become perforated and adapt the toroidal shape (see figure 12b, second from left). From 
a comparison between the experimental image and the simulated 3D structure we observe that cross section of the B 
block in the centre of the pore (figure 12b, utter right) looks very similar to figures 1C and D in Shin et al [16]. Since 
the structure for larger R and the same parameters is L\\ we conclude that this computational system gives a very 
close match. 



2. Other numerical studies 

The differences in methodology and choices of the parameter range complicate the comparison to results of other 
computational studies in general. Although most of these studies where not aimed at modelling experimental systems, 
we note that none of them have reproduced the experimental toroidal-like structure. [16] Moreover, the dynamics in 
our method is aimed at following the experimental pathway. Here, we focus on a qualitative comparison to the two 
most extensive MC studies for symmetric diblock copolymers in nanopores of three different radii: Feng et al [23] and 
Chen et al [22]. 

In Feng et al [23] xN is relatively low and stacked discs (L±_) arc found for neutral pores, as expected. The region 
where L±_ is stable shrinks for increasing pore radius. For large surface interactions, concentric cylinders (Lu) arc 
found. The only alternative structures for intermediate surface fields are PL (R/Lq « 0.8) and Lj_ u (R/L ~ 1.0). 
These results are fully in line with the findings in this article (see results section). Moreover, for R/L « 1.0 and an 
increased xN, also helices and double helices were obtained for varying surface field strengths. We also find these 
structures but for somewhat smaller R. Previously, we identified these structures as kinetically trapped or stabilized 
by clongational forces. The particular surface field strength and pore radius at which this effect plays an important 
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role may be expected to depend on the strength of the block-block interface, ie the value of \N. Moreover, also the 
slight asymmetry of our own block copolymer may lead to a shift of the parameter space where this effect is present. 
In Chen et al [22], the value of X-/V is considerably higher and the pore radii very small. For small to zero surface 
fields several structures are found to be stable: stacked discs L±, helices H and perforated (concentric) cylinders PL. 
Strangely enough tilted stacked discs are absent, which may have to do with the absence of extensional stresses for 
the considered pore length. Previously (previous section) we found that for small radii and our choice of %./V that the 
free energy difference between the helical and perforated lamellar structures can be small. Our results also suggested 
that the helical structure can be formed via the merging of perforations in the perforated structure, and subsequently 
that PL is an intermediate structure in the D — > PL — > H evolution. The fact that here PL is found for small 
surface fields is therefore not completely unexpected. For large surface fields, L| is stable. The only true discrepancy 
between this study and ours is the stacked circle structure found for D = 26 [22] close to the stability region of the 
L|| structure. This structure, similar to Lu but with short bounded instead of an infinite cylinder in the pore centre, 
could originate from packing problems in combination with a relatively large surface field. Although the bulk distance 
Lo was not determined explicitly, we can take Lq — 16 from counting the number of stacked discs in the cylinder, 
and find R/Lq = 0.81. We refrain from further speculations about the origin of this structure, and conclude that it 
is either an exotic structure or due to the higher \N . 

Finally, we briefly consider the results of a nanopore study that employed dissipative particle dynamics (DPD). In 
Feng et al [42] a A5B5 block copolymer system is simulated, both in neutral pores and pores with surface fields. 
The pore wall is not included as a geometrical constraint, but as soft core DPD-beads with strong repulsion. The 
effect of extensional forces is not considered. The authors find that L±_ always forms in neutral pores, except for the 
smallest pore radius considered (R = 7 [42]) where a double helix DH was found. For pores with strong surface fields 
only L|| forms, except for the smallest radius where L± is obtained. For surface field strength between these two 
extrema, several structures are obtained and hopping between different states of the system is observed. Although 
not mentioned as such, the structures in Figure 5 [42] show that they can be arrested in alternative morphologies. 
In particular, helical structures were found to serve as intermediates. Although the parameter range in this study is 
different and limited, it confirms our findings that the system can rather easily adapt to incommensurablity in larger 
nanopores. Moreover, for the intermediate parameter range (small, nonzero surface fields) there is a delicate balance 
between entropic and energetic contributions to the free energy, and structural transitions are mediated by helicity. 
In our more systematic study we show that structures can sometimes be arrested in these transition (helical) states. 
Moreover, we showed that for very small pore radii perpendicular structures can be found for relatively large surface 
fields. 



IV. CONCLUSIONS 

We have used a DDFT method to study the behavior of a slightly asymmetric diblock copolymer S10-B12 confined in 
a cylindrical nanopore. This block copolymer forms lamellae in the bulk, and was chosen to model the experimental 
system in Xiang et al [17]. We find many structures that were also identified in nanopore studies for completely 
symmetric diblock copolymers, but also a few new structures. However, almost all exotic structures are metastable. 
Our aim here is to identify general rules for structure formation in nanopores. 

We have focusscd on the intermediate surface strength in the parameter space, where the effective surface field, 
determined by the difference of the two block-surface interactions £, is relatively small. Outside this region, curved 
concentric cylinders (in analogy with the thin film notation denoted by L\\) are stable based on energetic arguments. 
The sign of the surface field determines the block that is preferentially found to wet the pore surface: for positive £ 
the B block is found close to the pore surface, for negative £ the S block is preferential close to this surface. Moreover, 
for a small and vanishing surface field the elastic chain deformations present in the concentric cylinders are avoided 
by the formation of a stacked disc structure (denoted as L±). The transition point from parallel structures (L||) 
to perpendicular structures (L±) is located somewhere inside this region, and we find many alternative structures. 
Analysis of the free energy and the determination of a phase diagram of equilibrium structures shows that all of these 
structures except Lu, PL and L± are metastable, and we have considered possible causes in detail. For stronger 
surface fields, that in general favour surface reconstructions (parallel structures), commcnsurability of the pore radius 
R and the natural lamellar domain distance Lq is important to avoid large compression or extension of the chains in 
parallel structures. For weaker surface fields, favouring perpendicular structures, commensurability of the pore length 
L x and the natural lamellar domain distance L is important to avoid large extensional forces. The commensurability 
issues, together with sometimes extremely slow transition dynamics, were found to cause the formation of alternative 
structures. 

For stronger surface fields, phase separation is initiated by the presence of the interacting pore surface. As a result, the 
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pathway from the homogeneous melt to more stable alternative structures includes a first stage of parallel structures 
or surface reconstructions, even when L\\ is not equilibrium. Subsequently, when the driving force is relatively large, 
holes are formed (for larger R, due to undulations that form connections between different layers) and connect: a 
helical structure is formed at a larger scale. The most probable origin of this helicity is that helical structures can 
more easily accommodate structural reformations at a higher level by adapting their spacing and pitch than defected 
L± structures. At the final stages, the helical structures line up by defect movement: the formation and destruction 
of three-folded connections. We find that these helical structures cost few energy to the system. Along this pathway, 
the structure can get kinetically trapped depending on the frustration that it experiences. We find metastable Lu, 
structures with holes (denoted by PL, and similar to the mesh or catenoid structures in [23] and [22], respectively), 
helices or double helices for decreasing surface field strength. In general, the transformation of the intermediate 
structures to other structures can be slowed down, or even be halted, by the presence of extensional forces due to 
incommensurability along the pore. 

For weaker surface fields the influence of the confinement on the initial stages is less substantial. As a result the 
initial structure is very interconnected and defected. The mechanism for the removal of these defects, however, is very 
similar to the one described above: small helical domains and three-fold connections play an important role. Here, 
the extensional forces play a major role, and determines whether stacked discs, tilted stacked discs or helical-like 
structures (which are actually clusters of tilted stacked discs under opposite angles, connected via helical domains) 
are formed. 

We have also considered the experimental system of Shin et al [16], in which a new toroidal- like structure was observed. 
Our results suggest that for weaker surface fields and small pore radius the structures are very sensitive to pore radius 
variations. Moreover, these pore radius variations give rise to interfaces between the stable structures for constant pore 
radii that exist in parts of the pore. These interfaces act as the driving force for the transition to other structures. 
For a certain surface field strength (£ = 0.4) and constant pore radius, we find a cage-like structure for a pore 
radius close to the experimental value, and perforated concentric cylinders (that appear as concentric cylinders in the 
experimental imaging procedure) for large pore radii. Moreover, by including the roughness of the experimental pore 
into our calculation in the form of a step profile, we were able to obtain a 3D structure that is very similar to the 
experimental structure. 
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VI. CAPTIONS 



Table I: Radial distance between 5-Bor block-surface interfaces (denoted by D t , J2i A = ^ counted from the 
pore surface) in parallel structures for £ = —1.0 and £ = 0.7. The position R of the interfaces is determined by the 
condition 6 S {R) = 6 B (R). 

Figure 1: Structural diagram with existing simulation results. Different symbols denote systems in different seg- 
regation regimes: triangles up - a system with Nx ~ 15 and triangles down the same system with Nx ~ 50, 
[23] circles - a A%B% system with Nx = 16, [13] squares - a AiqBiq system with Nx ~ 100. [22] Open symbols 
represent stacked discs, closed symbols concentric cylinders and red symbols alternative structures. We have scaled 
the axes to make a comparison: the vertical axis shows the reduced pore radius R/Lq, with R the pore diameter and 
Lo the particular lamellar spacing in bulk, the horizontal axis |e|, the energetic interactions of A with the surface in 
the MC method. We have used the relation \ ~ 5(e/kT) [33] for the conversion of our energetic FH parameter. 

Figure 2: Nanopore structure diagram for varying pore radius (in R/Lq) and surface field (in kJ/mol). Stable 
structures for each R/Lq (the determination of which is based on the free energies for constant pore length L x , see in 
Figure 5) are denoted by a grey backgound. 

Figure 3: Representative 3-D nanopore structures for: a) Lu, b) Lj_, c) L^, d) L±\\, e) PL, f) H, g) DH, h) 
GH. The corresponding locations in the structure diagram are: a) (1.74,-0.3) and (0.93,0.7), b) (1.74, 0.1) and 
(0.70,0.2), c) (0.93,0.1) and (1.05,0.1), d) (0.93, -0.1), e) (0.70,-1.0) and (0.70,0.7), f) (0.70, -0.1) and (0.70,-0.2), g) 
(0.81,-0.5), and h) (0.93,0.3), (1.05,0.5) and (1.28,0.4). Isodensity surfaces of the S component for mean value are 
shown. 

Figure 4: Representative 3-D nanopore structures for: a) dL\\, b) dL±, c) dL^ 11 , d) L 1 ^ 1 /L 1 ^ 1 , e) f) dPL, g) 

dDH, h) dGH . The corresponding locations in the structure diagram are: a) (0.93,-0.5), b) (2.21, 0.3) and (1.74,0.4), 
c) (2.21,-0.1) and (1.05,-0.1), d) (1.28, 0.1), e) (0.93,0.0), f) (1.05,-0.3) and (1.74,-0.2), g) (0.81,-0.1) and (0.81,-0.3), 
and h) (1.05,0.3) and (1.05,0.4). Isodensity surfaces of the S component for mean value are shown. 

Figure 5: Free energy (vertical axis) versus surface field strength £ (horizontal axis) for the reduced radii consid- 
ered: a) 0.58, b) 0.70, c) 0.81, d) 0.93, e) 1.05, f) 1.28, g) 1.74 and h) 2.21. The free energies associated with the 
structures in Figs 2 are denoted by open circles (for L x = 32) and open squares (for L x = 34 or 36). Similar structures 
are connected by straight lines to guide the eye. The open squares in a)-f) are associated with extra calculations 
(in all cases L±) for different L x : a) 34, b)-e) 36 and f) 34. In c) the free energy associated with the DH structure 
(£ = —0.5) is calculated for other £ (open diamonds). 

Figure 6: Nanopore phase diagram for varying pore radius (in R/Lq) and surface field (in kJ/mol). Stable structures 
for each R/Lq were determined using all free energy data in Figure 5, irrespective of the pore lengths L x . The absence 
of L|| for £ < (R — 0.58,0.70) and £ > (R = 0.70) does not allow for the analysis of the boundaries of the PL 
region. 

Figure 7: Kinetics of nanopore structure formation for selected pore radii and surface field strength (in kJ/mol): 
a) R/L = 0.70, £ = -0.4, b) R/L = 0.70, £ = 0.1, c) R/L = 0.81, £ = -0.5, d) R/L = 0.93, £ = -0.1, _e) 
R/L = 1.28, £ = 0.2, f) R/Lq = 1.74, £ = 0.3. The images show isosurfaces for the S'-block, for isosurface value 6s- 
The numbers indicate the dimensionless time steps TMS. 

Figure 8: The effect of extensional stress. Kinetics of nanopore structure formation for R/Lq = 0.70, £ = —0.5 
and varying pore length L x . a) L x = 32 (stable structure: H) and b) L x = 36 (stable structure: L±). The images 
show isosurfaces for the S-block, for isosurface value 6s- The numbers indicate the dimensionless time steps TMS. 

Figure 9: The effect of extensional stress for the double helical structure (R/Lq = 0.81, £ = —0.5). Upon a variation 
of the pore length L x a number of metastable structure clas_ses (denoted by DHi (i = 1 — 3)) were identified. The 
images show isosurfaces for the S'-block, for isosurface value 6s- Numbers indicate the pore length L x (in grid units). 

Figure 10: Evolution of the free energy associated with the structures in Figure 9. Numbers indicate the pore 
length L x . 
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Figure 11: Nanopore structures for 5*10512 (x = 0-6, L x = 32), matching the experimental system of Shin et al 
[16], a) for pore radii R = 19 (L x = 36), b) R = 9 (L x = 32) and surface field -1.0 < £ < 0.7 (in kJ/mol). Parallel 
structures for surface fields higher than the threshold surface field strength are not shown. The structures in c) are for 
constant surface field (£ = 0.4) and pore radius varying between 5 and 11 (L x = 32). For each radius, R, structures 
from two different viewing angles are shown. The images show isosurfaces for the S'-block, for isosurface value 6s- 

Figure 12: a) High isodensity image of the cage-like structure for R = 9 and £ = 0.4. b) Structures for the sys- 
tem of Figure 9 in a rough pore and £ = 0.4. The pore radius varies between R = 10 (center region) and R = 9 (outer 
regions). From left to right: isosurface for the S'-block (isosurface value 9s) after 1000 TMS, the same for 20000 TMS, 
isosurface of the B block (isosurface value 9b) after 20000 TMS, crop of the isosurface at 20000 TMS combined with 
an orthogonal slice of the B field in the centre of the pore. 
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